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We have studied experimentally the magnetic field-induced 
superconductor-insulator quantum phase transition in one- 
dimensional arrays ol small Josephson junctions. The zero 
bias resistance was found to display a drastic change upon 
application of a small magnetic field; this result was analyzed 
in context of the superfiuid-insulator transition in one dimen- 
sion. A scaling analysis suggests a power law dependence of 
the correlation length instead of an exponential one. The dy- 
namical exponents z were determined to be close to 1, and 
the correlation length critical exponents were also found to 
be about 0.3 and 0.6 in the two groups of measured samples. 

PACS numbers: 73.40.Gk, 73.23.Hk, 74.50.-Fr 



A ID array of small Josephson junctions (JJ) pro- 
vides an ideal testing ground for the T=0 quantum phase 
transition ||l|,||. Theoretically, a d-dimensional J J ar- 
ray can be mapped to the classical ((i-l-l)-dimensional 
XY model, however, the types of governing phase tran- 
sitions may not be the same in different dimensions 
In the 2D XY model, for instance, the ex- 
ponentially dependent correlation length of Kosterlitz- 
Thouless-Berezinskii(KTB) transition should lead to 
scaling properties different from those in a 3D model [|. 
While scaling properties in 2D J J arrays have been shown 
to exhibit an correlation length which corresponds to an 
underlying power law [^j^ , no ID arrays experiment has 
been reported to date. In this study we investigate the 
scaling behavior in ID JJ arrays, and report a power law 
dependent correlation length. In addition, correlation 
length exponents u are determined at approximately 0.3 
and 0.6. 

As shown in Fig. [^(a), the measured arrays were com- 
posed of Tshaped aluminum islands, whose sizes as de- 
fined by electron beam lithography were on the order of 
1 ixm. Each Al island has two tunnel junctions connected 
in parallel to its nearest neighbors, and forms a SQUID 
which is referred to as a unit cell. Being fabricated on 
the same chip, each group of measured arrays (denoted as 
group A and B) having different cell numbers, N , should 
have almost similar controlled parameters. Thus, length 
dependence of the phase transition can be investigated. 
The normal state resistances Rt of each cell for the ar- 
rays, as listed in Table 1, are determined at high bias, 
V > N{2A^ /e). The resistances for arrays in a group 
are quite similar, confirming the uniformity of the fabri- 



cated arrays. 

When magnetic field B is applied perpendicularly to 
SQUID loops with area A, the Josephson coupling energy 
Ej can be tuned periodically as Ej = Ej cos{7rAB /^o). 
The zero field Josephson coupling Ej is determined us- 
ing the r = Ambegaokar-Baratoff formula Ej = 
{A^ /2){Rq / Rt) with a superconducting gap A° of about 
200/ieV, and a resistance quantum Rq of h/Ae^ « 6.5kri. 
Accordingly, the Ej values are approximately 130//eV 
for arrays in group A, and 650/ieV for arrays in group 
B. From the SEM image one can estimate the junction 
area, and infer a junction capacitance C of 3.0 ±0.8fF 
by using a specific capacitance of 45 fF/ fim^ . With 
this capacitance, the charging energy for Cooper pairs, 
Ecp = (2e)V2C, is about 106±35/ueV, giving Ef/Ecp 
values of about 6 and 1.3 for arrays in groups A and B, 
respectively. 

The transport measurements were conducted in 
a dilution refrigerator at temperatures ranging be- 
tween 40mK and 1.5K. Samples were mounted inside 
an electromagnetic-field-sealed compartment within an 
oxygen-free copper holder. Any external noise was 
carefully prevented from reaching the sample compart- 
ment. The signal leads were filtered by 7r-section low- 
pass feedthrough filters, which were on the top of the 
cryostat and also with ThermoCoax cables from room 
temperature to mK temperature of sample compartment. 
To minimize the line frequency noise, a battery-powered 
preamplifier mounted on the top of cryostat was used. 
Furthermore, taking advantage of the common mode 
noise rejection, the entire measurement circuit was placed 
symmetrically with respect to its ground. Zero bias re- 
sistance Rq was taken from the slope of I{V) character- 
istics at a very small bias, and it was further confirmed 
by using a lock-in technique at a frequency of l.THz with 
an excitation of about 20nV. Sweeping a wide range of 
magnetic fields at 40mK, we determined the periodicity 
of magnctoresistance oscillation to be 9.15 Gauss. With 
this period AB, we denote the field as a dimensionless 
filling number / — B/ AB, which represents an average 
number of flux quantum in one cell. At integer values 
of /, the arrays are most conductive with i?o at a min- 
imum, while at half integer /-values, the arrays become 
most resistive with i?o at a maximum. 

Figure |l|(b) shows I{V) characteristics for array Al 
measured at / = ~0.5. For array Al at / = {i.e. 
the most superconducting curve), deviations of the su- 
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percurrent of consisting junctions are quite small, reaf- 
firming the uniformity of these arrays. However, even 
at /=0, the array is not truly superconducting, but has 
a finite zero bias resistance of about 0.9 kil. The su- 
percurrent decreases with / and diminishes at / ~0.45, 
above which the supercurrent-type structure turns into 
a Coulomb-blockade-type structure, with zero bias re- 
sistance reaching a maximum value of about 17 MQ at 
f—Q.5. The evolution from one structure to the other 
is best represented by Fig. |l|(c), which shows a smooth 
crossover from dip to hump structure with differential re- 
sistance Rd = dV/dl as a function of bias voltage, and 
with a flat Rd(V) curve separating the two limits. 

Due to a smaller Ej/Ecp value, the I{V) character- 
istics for array Bl, depicted in Fig. ||(a), have a higher 
Ro of about 100 kn at /=0 and 50 Gfl at /=0.5. The 
evolution of I{V) characteristics, from the supercurrent- 
type structure to Coulomb-blockade-type structure, for 
this array is quite different than that of array Al: As 
/ increases beyond /=0.20, a small Coulomb gap in the 
I{V) characteristic appears in the begining, signifying 
competition between Josephson coupling and Coulomb 
blockade of Cooper pair tunneling. This behavior can be 
clearly seen in Fig. ||(b), which shows, in addition to the 
simple dip and hump structures (as in Fig. |l|(c)), a coexis- 
tence of the dip and hump structures in the Rd{V) curves. 
From these Rd{V) curves, supercurrent and Coulomb 
blockade thresholds are plotted as a function of the fill- 
ing numbers and are shown in Fig. ||(c). Note that the 
two curves cross each other, a feature different from that 
reported in Rcf. ||l| with a considerably larger Ej/Ecp 
value (about 6.1) than those of our arrays in group B. 

The temperature dependences of zero bias resistance 
at various filling numbers for the arrays A2 and Bl are 
depicted in Fig. ^ For the most insulating case, the con- 
ductance fits standard Arrhenius form between IK and 
150mK with an energy barrier of about 120/ieV, which 
is very close to Ecp- This suggests that the dynamics is 
dominated by simple thermal activation of Cooper pairs, 
because the strength of the Josephson coupling is sup- 
pressed to the minimum. At even lower temperatures, 
shorter arrays show a saturation of resistance probably 
due to the finite size effect. For the most superconducting 
case {i.e. f—0), the resistance decreases with decreasing 
temperature and levels off at T <0.5K. At low temper- 
atures with / 7^0, Rq{T) increases even with decreas- 
ing temperature. A close inspection reveals a separation 
point Ro{T ^ 0) ~ Rq, above which Ro{T) curves go 
upward with decreasing temperature (leveling off for the 
case of shorter arrays) and below which Ro{T) curves 
simply level off. Similar results are reported by Haviland 
et al. [Q . We emphasize that this leveling off behavior 
is not due to any unwanted high-frequency noises in the 
measurement system for these reasons: 1. the leveling- 
off temperatures for the S and I sides are not the same, 
and 2. the leveling-off temperatures for the arrays in the 



two groups are not the same. This leveling-off behav- 
ior cannot be accounted for by heating effects because 
the measured resistance is almost the same for both AC 
and DC measurements. There are possible reasons that 
may bring in a finite leveling-off zero bias resistance: the 
finite size effect, the vortex macroscopic quantum tun- 
neling (MQT) 1^, and the charge Coulomb interaction 
]l0[ . The finite size effect can be ruled out because the 
leveling-off temperature does not change in arrays with 
different lengths. Even if it does involve, it should set 
in at the same temperature for both S and I sides, con- 
trary to the Rq{T) curves. The vortex MQT, observed 
in wire experiments has a pronounced effect in ID 
systems and can result in broadening of superconductor 
transition. However, MQT should be exponentially sup- 
pressed by dissipation, which is quantified as a ratio be- 
tween the quantum resistance and the junction tunneling 
resistance, Rq/Rt. Arrays in group A have a dissipation 
about five times larger than that of arrays in group B, 
and should have a diminishing vortex MQT. The fact 
that a higher leveling off temperature for arrays in group 
A suggests that MQT is not responsible for this behavior 
at this temperature. 

The last effect seems to be a reasonable one if the 
charge soliton picture is taken into account. While 
the origin model |Q describing phase transitions in ID 
Josephson junction arrays considered only an "on-site" 
Coulomb interaction, in reality, the long range interac- 
tion due to a nonzero island-to-ground capacitance Co 
cannot be ignored. For finite C/Cq, it has been shown 
PI that the charge will spread out and extend to a char- 
acteristic charge soliton length of A = y/ C/Cq. The 
result of the long range Coulomb interaction may lead 
to finite array resistances at low temperatures, driving 
even the superconducting phase to a metallic one [p^ . 
To investigate the effect of this charge soliton picture, 
we made nearby grounded electrodes along each array to 
control Cq. The distance between electrodes and arrays 
for samples in group A is 8 times longer than that for ar- 
rays in group B, leading to a smaller Cq and thus a longer 
soliton length. This is qualitatively consistent with our 
experimental finding that arrays in group A have higher 
leveling off temperatures than that of group B. 

For T > 0.5K, the sign of dRo/dT changes from pos- 
itive to negative value depending upon increasing filling 
numbers, signifying that the system undergoes a quan- 
tum superconductor-insulator transition. According to 
the theory of superfluid-insulator transition in ID sys- 
tems |lj], right at critical point / = /*, the resistance 
is linearly proportional to the temperature. Experimen- 
tally, /* is identified as the filling number where the ex- 
trapolation of Rq{T) curve passes Rq = 0, T = point. 
As a supplementary clue, we notice that at the base tem- 
perature, this critical filling number corresponds to an 
onset of Coulomb blockade threshold voltage (see Fig. 
||(c)). We thus interpret these phenomena as evidence of 



2 



a magnetic field-tuned SI phase transition with a critical 
filling number /* in our measured arrays. 

In a non-interacting model, the Hamiltonian of a ID 
array of small JJ can be mapped to a classical 2D XY 
model. Theoretically, the dimensionless coupling con- 
stant K, playing the role of the temperature in the clas- 
sical model, is related to the charging energy Ecp and 
Josephson coupling energy Ej as, K = Ej /2Ecp in 
the quantum system g . The 2D XY model has a KTB 
type transition. Below the transition temperature Tktb, 
the spins form vortex-antivortex pairs, while above Tktb , 
the pairs dissociate and the whole system becomes a vor- 
tex plasma Note that the topological spin vortex in 
the 2D XY model represents the phase slip in ID JJ 
arrays. In our system, the corresponding KTB transi- 
tion takes place at critical coupling energy E*j, which is 
achieved by the tuning of external magnetic field. In the 
region Ej < Ej, corresponding to the pair-dissociation 
phase at r > Tktb, the long range order in phase van- 
ishes, and phase fluctuations render insulating ID JJ ar- 
rays. According to the model, the transition takes place 
at E*j/Ecp = S/tt^ ~ 0.81 (s); however the experimental 
values are slightly larger than the theoretical one. De- 
spite the scattering of Ej/Ecp values in group A, the 
values in group B are well consistent to one other. The 
difference between the two groups can be explained with 
the effect of dissipation. At different strength of dissi- 
pation, it is expected that the critical point of SI phase 
transition will also be different ]l^ . 

The finite-temperature scaling law of quantum 
phase transitions asserts that 0{k,uj,K,LT) — 

L'^°/''0{kL](\ujLT,LT/£.T), where Lt = hP is a finite 
length on the imaginary time axis. Some of the terms can 
be neglected: the wave vector k is assumed zero, the scal- 
ing dimension do = 2 — d = 1 ||l^ , and w = in DC mea- 
surements. We thus obtain a concise finite-size scaling 
form for zero bias resistance, i?o(/, T) = T^^'-Ro{l/T^''), 
where ^ is a function of /. To examine the correla- 
tion length dependence, we rewrite the scaling form as 
RoiS,T) = T^''-Ro{Ti/T) where 5 ^ {K ~ K*)/K* - 
if ~ f*) is the distance from the transition point and Ti 
are field dependent scaling parameters (see for de- 
tails). In this way, one obtains a clear power law depen- 
dence of Ti on i5 over one order of magnitude, suggesting 
a power law dependence of the correlation length. This is 
in contradiction to the exponential dependence predicted 
in KTB transition and other theories JT^ . 

Based on this analysis, the correlation length ^ can 
be written as ^ = \5\~'^ where v is a critical exponent. 
Finally, tjhe scaling law can be written with a scaling 
function Rq as 

R,{f,T)^T'/^-xRo(^-^y (1) 

The exponent vz is determined from the slope of the 
power law fitting. For the seven arrays measured, we 



find that the uz are different varying from 0.30 to 0.60 
for different arrays. Nevertheless, the same values of lyz 
are obtained from both S and I sides. Assuming a z value 
of unity and using Eq. (|l]) with i/z and /* obtained from 
above, the scaling curves shown in Fig. ^ can be obtained. 
In addition, the scaling function form is found to be 
Ro{x) oc e"^, with a ~ ±10 and ±2.5 (in unit of Kl/''^ 
'+' for I-phase and ' — ' for S-phase), for arrays in groups 
A and B, respectively. This form for resistance is similar 
to results deduced from variable-range hopping(VRII) for 
the Bose glass |l^] . However, the VRH mechanism should 
not be accounted for our system since the exponent on 
T is larger than 1. The fact that the scaling function 
in one phase is symmetrical to that in the other phases 
suggests that the S phase and the I phase plays a dual 
role at zero bias. To refine the critical filling number 
and the scaling exponents, one begins with a trial scaling 
form, namely, i?o(/,T) = AT^^ exp(K(/)/Ti/''^). By 
noting that A is /-independent constant and that k is 
zero at /*, the vz and /* values can be unambiguously 
determined. The only adjustable parameter, z, can then 
be determined from the scaling curves. This parameter 
can be determined to a reasonable accuracy; a smaller z 
would give better scaling on the S-side, whereas a larger 
z would result in better scaling on the I-side. This fine- 
tuning procedure slightly modifies vz and /* values, and 
gives a z value of about 0.85±0.05 for array Bl. 

To summarize, we observed a magnetic field-tuned SI 
phase transition in the system of ID small Josephson 
junction arrays, and have, for the first time, made scal- 
ing analysis on such a system. Near the critical point, the 
Ro{T) scaling analysis indicates a power law dependent 
correlation length. The exponents vz are found to be 0.3 
to 0.6 with z close to 1, implying an isotropy in spatial 
and time dimensions. The value of correlation length ex- 
ponent ly contradicts what expected under current theory 
of ID Boson system, v > 2/d and i/ = oo jl^. The re- 
gions of lower temperature, where probably undermined 
complicated physics to make scaling fail, is away from the 
scope of the non-interacting model. These results suggest 
that certain important physics has not been unearthed in 
this system. 
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FIG. 3. The Ro{T) as a function of the fiUing number / for 
the array A2(a) and Bl(b). At about 0.5K above, the array 
displays an /-tuned SI transition, whereas at low tempera- 
tures the Ro{T) curves in the S-side level off or rise. 



FIG. 4. Ro{f,T) Scahng curve for array Bl for 
0.5K< T <1K and 0< / <0.5. The inset shows the scal- 
ing curve for array A2. 



TABLE I. Some important parameters of the measured ar- 
rays. The arrays that are fabricated on the same chip and 
thus have similar junction parameters, are categorized into 
one group. Note the closeness between the critical filling num- 
bers /* and the correlation length exponents u for arrays in 
a group. 



Sample 


Al 


A2 


A3 


Bl 


B2 


B3 


B4 


N 


49 


29 


14 


100 


75 


50 


30 


Rt 


0.9 


1.1 


0.9 


5.0 


4.7 


4.3 


4.9 


r 


0.41 


0.41 


0.45 


0.20 


0.25 


0.26 


0.27 


vz 


0.3 


0.3 


0.3 


0.6 


0.6 


0.5 


0.5 


E'j/Ec 


2.0 


1.6 


1.1 


1.05 


0.97 


1.02 


0.88 



FIG. 1. (a)The SEM image of an ID JJ array. The over- 
lapping areas between the T' shaped islands are the tunnel 
junctions. The scale bar at the bottom of the image is l^m. 
(b) Evolution of IV characteristics for arrays A2, from super- 
conducting behavior to insulating behavior, at selected filling 
numbers between 0.0 and 0.5. (c) Dynamic resistance Rd as a 
function of bias voltage show a crossover from superconduct- 
ing (bottom, /=0.42) to insulating behavior (top, /=0.50) 



FIG. 2. IV characteristics (a) and differential resistance 
Rd{V) (b) for array Bl for / = ~ 0.5 at 40mK. Notice 
the coexistence of the hump and dip structures; this is not 
seen in arrays in group A (c./. Fig. 1(c)). The curves in (b) 
are shifted for clarity, (c) The /-dependence of supercurrent 
(solid square) and threshold voltage (open diamond) for array 
Bl. Both are converted to the energy scale, i.e. Tilcj'ie. and 
eV . The supercurrent is magnified 100 times. 
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